Anisotropy of the gap parameter in the hole-doped cuprates 
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The structure of the gap parameter (A^) for the hole-doped cuprates has been studied in the 
paper. In the first step, the equation for the anisotropic anomalous thermal average (t/ik) has 
been solved close to the Fermi energy. The obtained results indicate the (/ik-amplitude crossover 
from the d-wave BCS to non-BCS behavior, as the electron-electron-phonon potential increases. 
In consequence, the antinodal part of the gap parameter is very weak temperature dependent and 
above the critical temperature (Tc), it extends into the anomalous normal state to the pseudogap 
temperature. On the other hand, the values of At, which are close to the nodal part, are temperature 
dependent and they disappear for T > Tc. The presented results strongly suggested that the 
pseudogap in the hole-doped cuprates is connected with the antinodal part of 2 At. The model has 
been tested for the YBa2Cu307_5 superconductor. It has been shown that the theoretical results 
agree with the experimental data. 



The gap parameter (Ak) is strong anisotropic in the 
cuprates [l| , ■ In the case of the hole-doped supercon- 
ductors, the d-wave symmetry is dominant 0], Q, 01 • In 
particular, for the optimally doped YBa2Cu307_5, the 
experimental data suggests 83% c?-wave anisotropy Q. 

The issue of the wave symmetry for the electron-doped 
compounds is still unclear. Some experiments favor the s- 
wave symmetry 01 , others support the c?-wave anisotropy 
(penetration depth tricrystal photoemission 

El' fl2l- R aman scattering and point contact tun- 
nclin g Urn). Additionally, the point contact tunneling 
data [14|| and the penetration depth experiment [TBj sug- 
gest the energy gap crossover from the d- to s-wave sym- 
metry, as the doping changes from the undcrdoped to 
overdoped region. 

The gap parameter for the hole-doped cuprates can be 
most precisely measured by using the ARPES method 
[T6| . So far, the obtained results have been interpreted 
in the framework of the two different approaches. 

In the first case, the difference between the doping 
and temperature dependence of the gap in the nodal and 
antinodal region suggests that the pseudogap and the su- 
perconducting gap are independent 17 1, [131 . Additional 
support comes from the strong deviation from the stan- 
dard d-wave form of the energy gap in the underdoped 
region. The above fact is interpreted as composing of the 
d-wave superconducting gap and the remnant pseudogap 

In the second case, the pseudogap is considered as the 
precursor of the superconducting gap [2l|, [2^. The d- 
wave symmetry deviation in the underdoped region is 
connected with the existence of the high-harmonic pair- 
ing terms [2^ . 



In the presented paper, we have studied the anisotropy 
of the gap parameter for the hole-doped superconduc- 
tors in the framework of the recently introduced theory 
In the above model, the pairing 
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mechanism is based on the electron-phonon (EPH) and 
electron-electron-phonon (EEPH) interactions. It has 
been shown, that for strong EEPH coupling and T < Tc, 
the average gap parameter (Ajot) is very weak temper- 
ature dependent, and up to the critical temperature ex- 
tends into the anomalous normal state: (i) in the s-wave 
case to the Nernst temperature, (ii) for the d-wave sym- 
metry to the pseudogap temperature. We boldly under- 
line that the discussed model explains well the experi- 
mental dependence of the ratio i?i = 2A[^l/kBTc on 
doping for the reported superconductors in the terms of 
the few fundamental parameters calculated by using the 
ab initio methods [13]. The symbol a[^1 denotes the 
average gap parameter close to temperature of the zero 
Kelvin. 

In the presented approach, the effective Hamiltonian 
for the high-Tc superconductors takes the form [2^ : 



(1) 



The first term describes the non-interacting electrons: 



(2) 
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The band energy for the two-dimensional square lat- 
tice can be written as: Sk = — t7(k), where t de- 
notes the nearest-neighbor hopping integral and 7 (k) = 
2 [cos (kx) + cos (ky)]. The symbols cj^.^ and Ckcr are the 
creation and annihilation operators for the electron with 
momentum k and spin a. The electron-phonon and 
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electron-electron-phonon effective terms are given by: 



2jYp '^ki-o-'^-kicr'^-k2ffCk2-o-j 



kik2(T 



and 



EC^ki~k4 , If 
2 A^Q <T C-k3 a Ckg -IT '^-k4 



ki ~k4(T 



(4) 

The parameter A^q is the normahzation factor: iVo = 
1 / '^ the symbol X^k " denotes the sum over the states 
for which: \e\^\ < ujq. The quantity Wq denotes the char- 
acteristic phonon frequency (ujq is of the order of Debye 
frequency). The functions Vkjk2 ^nd ?7k4r^k4 represent 
the pairing potentials with the d-wave symmetry: 



and 



Fk,k2 =^'y(ki)r;(k2), (5) 



f/ki^k4 = Ur^ (ki) 77 (k2) V (kg) V (k4) , (6) 



where: ?y (k) = 2 [cos (fc^) — cos (fcy)]. The anomalous 
thermal average is defined as: ^ko- = (c-ko-Ck-o-)- 

The Green function for the superconducting state has 
been calculated using the operator (P) [g^- The result 
takes the form: 



((ckt|c-ki)),^„ = 



i^Ek:w^kk2'/'k2 

' Ml 



(7) 



where Wkik2 can be written as: Wkik2 = H:ik2 + 
l^Ek3k4 f^ki~k4'^k3 0k4' and (/«k = 0ki- The symbol: 
bJn = {2n — 1) denotes the n-th fermionic Matsubara 
frequency; /3 = l/fcsT, where fc^ is the Boltzmann con- 
stant. The quantity takes the form: 



(8) 



[-kf:^-^-r.}j (^E^kk20k2) 

_^+^(]^&K3<^^3)(^|:.k4c)_ • 



On the basis of Eq. (O, we have derived the thermo- 
dynamic equation for the anomalous thermal average: 



(i;;&(i^O0k.) 




FIG. 1: The amplitude of the anomalous thermal average 
close to the Fermi energy for selected values of the potential 
U. We have assumed: V = 0.02t and luq = 0.3t. 




FIG. 2: The phase of the anomalous thermal average close to 
the Fermi energy. We have assumed: V = 0.02t and coo = 
0.3t. 



where: 



tan 



Xk 



2iy^el + A/, 



(10) 



v + u 



X '7(k)xk, 



■2 (Pko 



]VoE^Kk3)0L 



(9) We notice that Eq. ^ has been obtained with help of 
the relation: ((^l^))ii^„ = (^-S)- In particular: 

^ E„ I32+;|tmJ = ^y;?^ tan y^TA^) ^ 

The Eq. © has been solved for 7000 points close to 
the Fermi energy. The numerical result can be obtained 
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FIG. 3: The amplitude of the anomalous thermal average close to the Fermi energy for selected values of U and the temperature. 
In particular, the figure (A) is for (7 = 0, (B) for U = 0.02t, (C) for U = 0.05t, and (D) for U = O.lt. We have assumed: 
V = 0.02t and cjo = 0.3t. 



with help of the self-consistent iterative (SCI) procedure 
originally used for solving the Eliashberg equations de- 
fined on the imaginary axis l2i,|3i,[3l|,[33,[3l,[33, 
m, Hi], [13 • We clearly underline that SCI procedure 
is very appropriate for the large sized system of the non- 
linear equations. The stability of the numerical solutions 
is ensured by replacing the initial values with the calcu- 
lated values and repeatedly running the computer pro- 
gram. The computational errors can be controlled by 
solving the set independently of each in the framework 
of the Fortran and C programming language. 

In Fig.m we have presented the amplitude of the 
anomalous thermal average (|0k|) close to the Fermi en- 
ergy for different values of the potential U . The tem- 
perature takes the very low value {ksT = 10~H). The 
obtained results show that beyond the nodal regions the 
low values of |(/)k| strongly increase with the growth of U . 
For high values of U , the observed transition from the d- 
wave BCS to non-BCS behavior results in coming into 
existence the nearly-constant value of the amplitude in 
wide surroundings of the antinodal regions. In the nodal 
regions appear the characteristic "icicles" created by the 
wide range of j^kl values. 

In Fig. m we have plotted the phase of the anomalous 
thermal average ((/Sk)- In contrast to the amplitude |(/)k|, 



the phase is independent from the parameter U . Addi- 
tionally, it has been shown that (/7k assumes two different 
values. 

The temperature dependence of |0k| for selected val- 
ues of U has been presented in Fig.[3l In the first step, 
the case of J7 = has been considered (Fig. [3] (A)). It is 
easy to see tliat the amplitude decreases with the grow 
of the temperature and disappearances at Tc — 0.0231i. 
In other cases ( Fig. [3] (B)-(D), where U ^ 0), the tem- 
perature dependence of |</)k| sharply differs from the pre- 
diction based on the pure d-wave model. In particular, 
for < T < Tc the antinodal regions of |0k| are very 
weak temperature dependent and above the critical tem- 
perature extend into the anomalous normal state to the 
pseudogap temperature. In contrast to the behavior of 
the antinodal values, the nodal regions of the amplitude 
strongly decrease with the grow of the temperature. 

Now, we can calculate the gap parameter: Ak = 



?7k \(l>v 



^ + c/|^k||<Ak| 



The form of the function Ak 

has been presented in Fig.[31 In particular, we have se- 
lected the gap parameter for J7 = 0, close to the zero 
temperature (fc^T = 10~H), and the gap in the range 
of high value of the clectron-electron-phonon potential 
{U ~ O.lt) for two following temperatures: fc^T = 10~^t 
and fcsT = 0.2339i. 
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FIG. 4: The gap parameter close to the Fermi energy for 
selected values of U and the temperature. In particular, the 
figure (A) is for f/ = 0, (B) and (C) for U = O.lt. We have 
assumed: V — 0.02t and ujo = 0.3t. The figures on the right 
show the general form of the gap; the blue surface represents 
the simple d-wave gap. The shape of the gap parameter in 
the nodal regions has been in detail presented in the figures 
on the left; the pink-blue surface represents the simple d-wave 
gap. 



In the case of J7 = 0, one can easily see that the gap 
parameter very well reconstruct the shape of the simple 
d-wave gap (A|^'''' = A?7(k)). However, the function At 
possesses enough complicated internal structure. We no- 
tice also that in the nodal region the gap closely adjoins 
to the A^J''' surface. Of course, the function Ak decreases 
in the typical way with the increase of the temperature 
and assumes the zero value at Tc — 0.0231t. 

For the high value of the clectron-clectron-phonon po- 
tential and fc^r = 10~^t, the amplitude of the gap strong 
increases and the internal structure of Ak becomes much 
straighter (see Fig.S] (B)). In the nodal region the gap 

agrees relatively well with the A^.''-' surface. 

In the case oi U = O.lt and ksT = 0.2339t, one can 
observe anomalous behavior of the gap (Fig. |4] (C)). 
In the antinodal region Ak has comparable values as for 
ksT = 10~'^t. From the physical point of view the ob- 
tained result indicates that the pseudogap exists in the 
antinodal region. On the other hand, in the nodal region 



the gap parameter forms the characteristic area of the 
very low values. 

The model has been tested on the example of the 
YBa2Cu307_5 compound [3^. It has been assumed that 
the nearest-neighbor hopping integral and the character- 
istic phonon frequ ency take the values 250 meV and 75 
meV, respectively [s^, [40l |. 

On the basis of the relations: Tq = Tc {V) and 
T* = T* {V,U), where Tc denotes the critical temper- 
ature and T* is the pseudogap temperature, the pairing 
potentials have been uniquely determined. Then, the 
equation ^ has been solved, and the values of 2A (a) 

have been calculated, wherein a = arctan • 

The theoretical and the experimental results have been 
compared with each other in Fig. [5] It can be easily seen 
that the presented model in the correct way predicts the 
course of the experimental data. 

In the paper, we have tested the microscopic model 
for the high-Tc superconductivity, which is based on 
the electron-phonon and electron-electron-phonon inter- 
action. In particular, the structure of the gap parameter 
has been studied. The presented results have been ob- 
tained using the thermodynamic equation for the anoma- 
lous thermal average. 

It has been stated that for U ^ 0, the gap parameter 
very well reconstruct the d-wave behavior. However, Ak 
possesses complicated internal structure. 

In the opposite case (the large value of U) the behav- 
ior of the gap parameter is anomalous. In particular, the 
antinodal part of Ak is very weak temperature depen- 
dent and above Tc it extends into the anomalous normal 
state to the pseudogap temperature. On the other hand, 
the nodal part is temperature dependent and for T > Tc 
the values of the gap disappear. The above results sug- 
gested that the pseudogap in the hole-doped cuprates is 
connected with the antinodal part of the energy gap. 

In the last paragraph, let us notice that the presented 
model has been tested for the YBa2Cu307_i compound. 
Good agreement between the theoretical predictions and 
the experimental results has been achieved. 
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FIG. 5: The dependence of the energy gap (2A) on the angle for the YBa2Cu307_4 compound. The circles represent the 
theoretical results. The squares are the experimental data [3q |. 
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